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Abstract 

The coalescence of liquid drops induces a higher level of complexity compared to the classi- 
cal studies about the aggregation of solid spheres. Yet, it is commonly believed that most findings 
on solid dispersions are directly applicable to liquid mixtures. Here, the state of the art in the 
evaluation of the flocculation rate of these two systems is reviewed. Special emphasis is made on 
the differences between suspensions and emulsions. In the case of suspensions, the stability ratio 
is commonly evaluated from the initial slope of the absorbance as a function of time under diffu- 
sive and reactive conditions. Puertas and de las Nieves (1997) developed a theoretical approach 
that allows the determination of the flocculation rate from the variation of the turbidity of a sam- 
ple as a function of time. Here, suitable modifications of the experimental procedure and the 
referred theoretical approach are implemented in order to calculate the values of the stability ra- 
tio and the flocculation rate corresponding to a dodecane-in-water nanoemulsion stabilized with 
sodium dodecyl sulfate. Four analytical expressions of the turbidity are tested, basically differing 
in the optical cross section of the aggregates formed. The first two models consider the processes 
of: a) aggregation (as described by Smoluchowski) and b) the instantaneous coalescence upon 
flocculation. The other two models account for the simultaneous occurrence of flocculation and 
coalescence. The latter reproduce the temporal variation of the turbidity in all cases studied (380 
< [NaCl] < 600 mM), providing a method of appraisal of the flocculation rate in nanoemulsions. 

Keywords: 

Aggregation, Coalescence, Flocculation rate, Nanoemulsion, Stability ratio 



Contents 

1 Introduction 

1.1 Flocculation of solid particles 

1.2 Flocculation of oil drops . . 



'Corresponding author 

Email address: krahn9ivic.gob.ve (Kareem Rahn-Chique) 
'Tel: +58-(212) 5041980. Fax: +58-(212) 5041915. 



Preprint submitted to Advances in Colloids and Interface Science 



March 4, 2013 



2 Theoretical background lfl 

2. 1 Theoretical models for the estimation of the turbidity curve 11 

2.1.1 Aggregation of emulsion drops (model Ml) 11 

2.1.2 Complete coalescence of the drops after contact (model M2) 13 

2.1.3 Partial coalescence of flocculating drops (model M3) 14 

2.1.4 Modification of Ml to account for a finite coalescence time (model M4) [14 

3 Experimental 15 

3.1 Materials 15 

3.2 Phase diagram [NaCl] vs. [SDS] 15 

3.3 Nanoemulsion preparation 16 

3.4 Nanoemulsion characterization 16 

3.5 Experimental determination of stability ratio ( W) 16 

3.6 Evaluation of W exp 17 

3.7 Evaluation of Wu u 17 

3.8 Fitting of the turbidity curve 17 

3.9 Polydispersity [17 

4 Results 

4. 1 Variation of turbidity as a function of time 

4.2 Fitting of the turbidity curve with theoretical models 

4.3 Effect of the initial polydispersity on the evaluation of kf 

5 Discussion |2^ 

6 Conclusions |23| 

7 Acknowledgements IH 

8 Appendices 24 

8.1 Appendix A: Alternative procedure for evaluation of Eqs.[l4][l5][T8] and[T9] 24 

8.2 Appendix B: Balance of aggregates for a flocculating-coalescing emulsion 25 

8.3 Appendix C: Calculation of errors 27 

8.4 Appendix D: Comparison between Eq. QUand Eq. |39](Model Ml) |27 

9 References 



1. Introduction 

The absolute coagulation rate is an important parameter used to characterize the coagula- 
tion kinetics of colloidal systems. Knowledge of the stability and flocculation kinetics of these 
systems is essential for different applications in chemical, food, pharmaceutical, and material 
industries, environmental engineering, biotechnology, and nanotechnology |Q], H 0] • 

Smoluchowski [5] considered the kinetics of irreversible flocculation of a homogeneous dis- 
persion, initially composed of single particles suspended in a liquid. The number concentration 
of aggregates of k primary particles (size k) existing in a dispersion at a given time («t(f)) results 
from a balance between the aggregates produced by the collisions of smaller clusters of sizes i 
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and j (i + j — k), and the number of aggregates of size k lost by a collision with aggregates of 
any other size: 
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The kernel of the Eq. Q]is the set of coagulation rates of aggregates of size i and j, Assuming 
equal fiocculation rates for all aggregates of different size (fey = kf), Smoluchowski deduced a 
simple analytical equation for the change in the total number of aggregates as a function of time 

(« a ): 

(2) 



(1 + kffiQt) 

Here, «o stands for the initial number of aggregates (n a (t = 0) = no)- The fiocculation rate sug- 
gested by Smoluchowski, kj = 4kgT/3ri (where: kg, T and rj stand for the Boltzmann constant, 
the absolute temperature and the viscosity of the external phase, respectively) solely considers 
the Brownian movement of the particles. 

Derjaguin-Landau-Verwey-Overbeek Hill HI, introduced the effect of the attractive (V A ) and 
the electrostatic (Ve) potentials in the process of fiocculation. According to their theory (DLVO), 
the total potential of interaction between the particles, ( Vt = Ve+Va) shows a primary minimum, 
a repulsive barrier, and a secondary minimum. Fiocculation in the primary minimum is assumed 
to be irreversible due to the strength of the van der Waals interaction at short distances. The onset 
of destabilization occurs when the maximum of the repulsive barrier of the potential is equal to 
zero (Vt = and j^- = 0, where ry is the distance between the particles i and j). Nonzero 
repulsive barriers provide kinetic stability to the system with respect to irreversible fiocculation 

To account for the effect of a repulsive barrier in the equation of Smoluchowski, Fuchs 111 ill 
introduced the concept of a stability ratio, W. Later on, McGown and Parfitt fl2ll modified it 



to include the attractive interaction in the calculation of fast aggregation, and Spielman 111311 
introduced a correction due to the hydrodynamic interaction: 
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where, is the radius of particle k. The pre-exponential factors of Eq. [3] account for the effect 
of the hydrodynamic interactions on the particles movement. 

1.1. Fiocculation of solid particles 

A comprehensive compendium of the early experimental techniques and theoretical approaches 



used to quantify the fiocculation rate of suspensions was published by Sonntag and Strenge [ 14] 
in 1987. A formal treatment of the hydrodynamic aspects of particle diffusion and electrokinetic 
effects can be found in the books of Russel et al. II 1 3f| . van de Ven 11611 . Lyklema 11711 . Hunter 



0181], Adamson [121, Elimelech et al. [20] and Dhont (21J]. Up to the 1980's most of the research 
on this subject dealt with: a) the testing of constant-charge and constant-potential solutions of 
the Poisson-Boltzman equation required for the calculation of the electrostatic interaction be- 
tween the particles, b) the appraisal of hydrodynamic interactions, c) the differences between the 
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electrostatic potential of the particles at the interface, and the zeta potential deduced from mo- 
bility measurements, d) the inclusion of non-DLVO forces between the particles, e) the effect of 
counterions on the stability of the suspension, and f) the consequences of reversible (secondary 
minimum) flocculation 11221 l23l 12411 . 

A review of the experimental methods currently employed for measuring aggregation rates 
was recently published by Gregory [25]. Lips et al. |26|, and also Lichtenbelt et al. I27I1 de- 
veloped experimental methodologies suitable for quantifying the rate of doublet formation (ku) 
with light scattering. The procedure reported in Ref. [26] allows obtaining the value of W from 
the initial slope of the turbidity as a function of time |26[ 28, 27, 29]. This methodology is very 
convenient since it only requires the use of a spectrophotometer. Reering and Overbeek deduced 
an analytical relationship between log W and log I s (where I s is the ionic strength of the solution) 
which was confirmed by the experiments ll30ll . Using numerical evaluations of the interaction 
potentials, Prieve and Ruckenstein reported a linear dependence between log W and the height 
of the repulsive barrier between the particles [Bill . 

In order to measure absolute coagulation rates, Young and Prieve rf29ll improved the small- 
angle scattering technique proposed by Lips and Willis [28] shortening the path length of light 
through the dispersion. For rapid flocculation occurring at high salt concentration, a stability 
ratio of 2.14 was obtained for a polystyrene latex of diameter d =486 nm. 

Sun, Liu and Xu J32, 33, 34]] made an extensive study regarding the necessary conditions for 
a reliable turbidity determination of the coagulation rate. The relative rate of turbidity change at 
short times: Rj = (l/To)(dT/di)o, is found to depend markedly on the wavelength of the incident 
light and the size parameter a = 27ia/A (where a is a characteristic length of the particle). If the 
extinction cross section of two singlets (2ct[) is less than that of one doublet (<x 2 ), then Rj > 0, 
otherwise Rj < 0. Depending on A, Rj can be equal to zero (blind point) if <x 2 ~ 2cr\ . In this 
case, the turbidity is insensitive to the aggregation process. In mixtures of particles with average 
radii R — 1 15, 266, and 500 nm, the blind point is found to depend on the size of the particles 
and on the fraction of each component i33ll . Furthermore, larger particles greatly mask smaller 
ones, inducing experimental errors. 

These authors 02l l33l 13411 demonstrated that the T-matrix method can calculate the optical 
cross sections of singlets and doublets accurately, and confirmed that the Rayleigh-Gans-Debye 
(RGD) 1 35] approximation is only reliable for very small particles: a < 3. When this condition 
is not fulfilled, the value of ku changes as a function of the wavelength. The value of A that 
maximizes Rt guarantees better measurements, but it is restricted by the additional constraint 
on a. Further uncertainties can be caused by unwanted contributions of the turbidity due to the 
forward scattering of light 113411 . 

Toward the end of the last century, the synthesis of functionalized latex particles of micron- 
size became widespread as well as the use of static (SLS) and dynamic (DLS) light scattering 
techniques. The latex particles have the size of a few hundred nanometers, and are almost neu- 
trally buoyant. As a result of these advances, the morphology of the aggregates was linked to 
the aggregation rate |36 , 37 , 3^, and some mathematical properties of the flocculation rate were 
quantified thanks to the dynamic scaling of the population balance equation for irreversible ag- 
gregation lf38l [39l l40l l4ll |42 . 43, 44]. From these studies, a general form of the flocculation rate 
between aggregates of sizes i and j was deduced J45ll : 
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Here K is the kernel for Brownian aggregation, B, ; is a geometric correction which accounts for 
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the fractal dimension of the aggregates, Wy is the stability ratio between particles i and j, and P (J 
is the variable that contains the details of the aggregation mechanism. 

Using single particle light scattering (SPLS), Holthoff et al. 114611 studied the time evolution 
of aggregates from doublets up to heptamers. For particles of d =580 and 683 nm, and I s = 1 
M KC1, the experimental data nicely fitted the predictions of Smoluchowski corresponding to a 



constant kernel, with &n = (7.3-7.1) xl0~ 18 m 3 s _1 [ 47] . Similar values were obtained using DLS 
and SLS techniques (6.9-6.1) xlCT 18 m 3 s~'. Moreover, a linear variation of log W vs. logI s was 
found |46i|. For a salt concentration of 250 mM, the stability ratio varied in the range of 2.5-2.8 



(11.9 < W < 14.2 for I s = 125 mM). Using DLS and SLS, Behrens et al. 048J] measured the 
absolute coagulation rate of doublet formation and compared it with the theoretical predictions of 
the standard DLVO theory. The dependence of the stability ratio with the pH, the surface charge 
of three particles (two carboxylic lattices (d =309, 104 nm) and hematite (d =122 nm)), and the 
ionic strength of the solution was studied. For the systems under consideration the predictions 
of the formulas corresponding to constant charge, constant potential and regulated charge do not 
differ much. It was concluded that for colloid particles of low surface charge (< 3 mC mr 2 ) and 
moderately low ionic strength (0.1 - 10 mM) the classical DLVO theory works quantitatively. 
However, when the ionic strength is increased to 100 mM, the maximum of the potential curve 
lies below 2.2 nm, and a deep secondary minimum develops. In this case, the slope of W vs. pH 
predicted by the calculations differs considerably from the one suggested by the experimental 
measurements. For these experiments, the value of W changed between 1 and 1000, and the fast 
flocculation rate for doublet formation (jd" st ) ranged from 2.4 xl0~ 18 to 3.5 Xl0 _18 m 3 ,s _1 . 

Nowadays, SLS and DLS techniques are commonly applied to study homo and hetero- 
aggregation in complex systems 1491 l50l Dill . These techniques allow obtaining the value of 
k\\ under reactive and diffusive conditions for solid particle dispersions. However, bigger ag- 
gregates are almost exclusively considered for topological studies, not for the calculation of the 
aggregation rate. 

1.2. Flocculation of oil drops 

Typical macro and nanoemulsions are kinetically stabilized dispersions of oil (o) in water 
(w). Unlike suspensions they are subject to coalescence, Ostwald ripening and phase inversion, 
besides flocculation and creaming. These phenomena occur simultaneously and are not inde- 
pendent of each other. Emulsion drops are usually of micrometer size and their coalescence 
dynamics includes the occurrence of deformation, capillary waves, and hole nucleation. The 
density of the oils varies between 0.6 and 0.9 gcm~ 3 , which generates a substantial amount of 
creaming upon aggregation and particle growth. For an overview of some of the problems related 



with emulsion stability, the reader is referred to the following books: [52, 53, 541 15511 . 

While big drops (R > 50pm) deform during flocculation, small drops (R < lpm) are reluctant 
to deform due to their high Laplace pressure 115 6L Whenever deformation occurs the draining of 



the intervening film between the drops takes a finite time I157U58L 15911 . Although the process of 
coalescence can be formidably complex [60], deformation allows distinguishing between floc- 
culation and coalescence. However, if two drops do not deform during their approach, the final 
state of the system depends on the interaction potential. A high repulsive barrier will promote 
secondary minimum flocculation as in the case of solid particles, but a small or non-existent 
barrier, will lead to a steady approach until coalescence occurs. In this case, it is not possible to 
determine when the process of flocculation finishes and the one of coalescence begins. Hence, an 
experimental quantification of a flocculation rate under diffusive conditions necessarily includes 
the process of coalescence 16111 . 
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The deformation of the drops depends on their size, hydrodynamic interactions, colloidal 
forces, and surface forces |57l 62]. According to recent emulsion stability simulations (ESS) an 
approximate threshold for drop deformation occurs around 2 .5 /J.m 

Drops with 

radii smaller than 2.5 fim will mostly behave as non-deformable particles. 

In the absence of surfactants, oil drops show a substantial electrostatic potential. This surface 
charge is commonly ascribed to the preferential adsorption of hydroxyl groups at the oil/water 



interface [65, 66| |67J]. Hence, it depends on the pH and on the total interfacial area of the 
emulsion. 

Using DLS, Sakai et al. 1 68] observed that the radii of benzene droplets dispersed in water 
increase in a geometrical series with a common ratio of (V/3V ) 1 ^ 3 , where V is the total volume 
of the emulsion (V = V w + V ), and V, is the volume of phase i. This result was explained 
supposing a multi-step-growth model which assumes that primary droplets first aggregate into 
a cluster, and then merge (coalesce) making a secondary droplet. Similarly, secondary droplets 
aggregate and then merge into a tertiary droplet, and so on. The rate of the process depends on 
the mean free path of the drops. 

Sakai et al. ll69ll studied the evolution of alkane-in-water nanoemulsions (R = 30 nm, <f> - 

2 xl0~ 4 , where <p is the volume fraction of oil) in the absence of surfactants. According to 
the theory of Lifshitz, Slyozov and Wagner (LSW theory) |70l 71], the Ostwald ripening rate 
(dR 3 jdt) is proportional to the solubility of the oil, and the interfacial tension of the drops. Our 
group carried out simulations on the behavior of the dodecane-in-water nanoemulsions studied 
by Sakai 17211 . In the absence of salt, the simulations reproduced the experimental behavior, 
showing that: a) the initial evolution of R 3 is only due to flocculation and coalescence, b) the 
slope of R 3 vs. t decreases with time approaching the LSW limit, and c) molecular diffusion only 
contributes to the augment in R after 800 s. At high salt concentration it was predicted that R 3 vs. 
t varies linearly as a product of flocculation and coalescence producing a polydisperse emulsion. 

Most of the works on emulsion stability concern macroemulsions. Using video enhanced 
microscopy, Holt et al. ll73l l74l 17511 studied the lifetime of a doublet (f^) in toluene in water 
emulsions (0 = 0.005-0.02, R ~ 2-3 nm) at low ionic strength (3-13 mM NaOH). At < 

3 mM there is no influence of I s over f^. At I s > 3 mM, the distribution of lifetimes splits 
into two parts: one at small contact times that do not depend on the value of I s , and another 
at large lifetimes which grow longer as the electrolyte concentration increases. It was found 
that the former lifetimes correspond to collisions between two particles in which the minimum 
distance of separation is longer that the position of the secondary minimum. The latter values 
of t c { correspond to secondary minimum flocculation. The depth of the minimum augments as 
the ionic strength rises, increasing the value of f</. Values of tj of the orders of several minutes 
were observed, as well as evidence of droplet coalescence inside doublets. For very high ionic 
strength (I s =0.1 M NaOH, W ~ 2), the decrease in the number of singlets closely follows the 
predictions of Smoluchowski. It was also observed that the drops tend to stick to the walls of the 
coalescence channel when I s > 0.1 M. 

Verbich et al. |[76ll studied the evolution of a 1 , 1 -dichlorodecane in water macroemulsion (<p = 
0.01, R ~ 2-4 fim, 10~ 4 < [NaOH] < 10~ 2 M). They found that the relative absorbance, Abs/Abs 
decreased with time during several days. The zeta potential of the drops also decreased unless 
SDS was added to the system. The value of the stability ratio was calculated using the time 
required for a given emulsion ([NaOH] = 10~ 3 M) to reach the same optical density of a fast 
flocculating one ([NaOH] =0.3 M) W ~ 48. The emulsions did not show a sensitive variation of 
coagulation rate with the ionic strength, which does not agree with the theory of Fuchs. Using 
optical microscopy it was observed that as the drops of 2 ymi disappear, the content of 6-10 
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fim drops increased gradually, and larger drops start to appear. Finally, the number of doublets 
consisting in one large droplet (12-16 jum) and one small droplet increased. Multiplets were 
observed after several days. 

Karaman et al. [77], Pashley 178], and Maeda et al. [79] suggested that dissolved gas has 
a key role in the stability of surfactant-free dodecane-in-water emulsions. Degassed oil sponta- 
neously emulsifies in water producing stable emulsions whose turbidity progressively decreases 
as a function of time. In apparent contradiction, the reintroduction of dissolved gases does not 
destabilize the emulsions already formed. The turbidity of emulsions not degassed was con- 
siderably lower, but also diminished with time. This was ascribed to the buoyancy removal of 
the large drops. Burnett et al. showed that the behavior of these emulsions is related to the 
freeze-pump-thaw treatment and not to the amount of gas dissolved [80, 81]. 

Ionic surfactants adsorb to the interface of the drops increasing their electrostatic charge 
1 661 82Tl . The charge depends on the equilibrium adsorption of the surfactant, which is a function 
of the surfactant concentration. At a certain concentration surfactants aggregate into micelles. 
This equilibrium is modified by the presence of salt in the aqueous media that decreases the 
critical micelle concentration and increases the surfactant surface excess. However, as the ionic 
strength increases, the additional charge is also partially screened by the salt counterions. 

High salt concentrations decrease the solubility of the surfactant in the aqueous phase. The 
surfactant may precipitate as hydrated crystals if is low lTs3h - In the most general case, high 
salt concentrations favor the formation of a separate surfactant-rich middle phase, or even the 
transfer of the surfactant to the oil phase in the form of inverse swollen micelles [84, 85, 8(J 871. 
The region of compositions for which a three-phase system exists at equilibrium is known as the 
balance zone. Within this region the interfacial tension produced by gemini surfactant attains 
ultralow values. Aveyard et al. [88] showed that a single chain surfactant like SDS also gener- 
ates a tension minimum in the balance zone, but it cannot reach ultralow values of the tension 
unless a significant amount of a cosurfactant is added to the system. Interestingly, the stability of 
emulsions in the balance zone is also minimum, changing in several orders of magnitude within 
a narrow region of compositions Hill 0|9l]]. This lack of stability has been the subject of 
several hypotheses. Recent works connect unstable emulsions with the spontaneous curvature of 
the surfactant at the interface, and the pronounced occurrence of Ostwald ripening. If the salt 
concentration is increased beyond the balance zone, the excess surfactant is transferred to the oil 
phase and the emulsion passes from o/w to w/o l84l[85l[86l.[83ll . 

Using silicon-in-water emulsions (0.34 < R < 3.3 pmi) stabilized with sodium dodecyl sulfate 
(SDS), Bibette et al. [92] demonstrated that a sharp destabilization threshold can be defined. In 
these experiments, dialysis bags and Dextran were used to set the value of the osmotic pressure. 
The threshold was found to be a function of the surfactant concentration, the osmotic pressure, 
and the droplet diameter. 

Schmitt et al. Ir93l l94ll studied the behavior of concentrated macroemulsions of alkane-in- 
water, noticing that Ostwald ripening dominates the evolution of the emulsion shortly after emul- 
sification. However, a few hours later, coalescence increases abruptly equalizing and surpassing 
the Ostwald ripening rate. The point at which the two rates of destabilization are equal, allows 
the quantification of the frequency of rupture of o/w/o films, which is a measure of the coales- 
cence frequency. This frequency is around 10 4 m~ 2 s~ l for the case of heptane-in-water emulsions 
stabilized with SDS. 

Aveyard et al. [95, 96] and Binks et al. [91 ] studied the flocculation transitions of dodecane- 
in-water and heptane-in-water macroemulsions (3 < R < 10 ymi) stabilized with non-ionic sur- 
factants. These emulsions are stable with respect to coalescence and Ostwald ripening, but un- 
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stable with respect to flocculation and creaming. The transitions between non-flocculated and 
flocculated states were induced by the addition of an ionic surfactant and salt. In this case, the 
surface charge of the drops depends on the partition of the ionic surfactant between the aque- 
ous solution and the interface. Using analytical forms for each potential contribution, the total 
interaction energy and the disjoining pressure between two deformable emulsion drops were cal- 
culated. This allowed the definition of regions of flocculation in plots of the ionic surfactant 
concentration vs. the total ionic strength. Additionally, the destabilization of the emulsion was 
quantified in terms of the creaming rate. The comparison between theory and experiment was 
very good. 

van den Tempel considered the process of coalescence of o/w systems as a first order re- 
action depending on the number of films between aggregated drops B97I1 . Using the equation 
of Smoluchowski and assuming linear aggregates, he was able to deduce an explicit analytical 
expression for the total number of aggregates in an emulsion, which depends on a flocculation 
and a coalescence rate. Borwankar et al. deduced a similar average expression for the number 
of aggregates, considering the coalescence of small floes 119 al . Fitting the theoretical expressions 
to the experimental data from soybean in water emulsions (<p = 0.30), flocculation rates of the 
order of 10~ 16 m 3 s~ l were obtained, along with coalescence rates between 7.4 xl0~ 6 and 2.5 

xio-V 1 . 

Reddy, Fogler and Melik 

GUM 

1 111 studied the simultaneous occurrence of flocculation 
and creaming by numerical evaluation of the respective population balance equation. This pro- 
cedure allowed explaining the cyclic shifts of the drop size distribution (DSD) toward higher and 
lower particles sizes, observed in paraffin/water emulsions (R ~ 1-4 pmi). Moreover, they re- 
produce the evolution of the DSD at different heights of the container. These distributions differ 
significantly from one another. 

A kinetic model for t he s imultaneous processes of flocculation and coalescence was for- 
warded by Danov et al. Ill 0211 . Most of their differential equations must be evaluated numer- 
ically. However, they demonstrated that if the coalescence rate of an emulsion is significantly 
faster than the flocculation rate, and big drops are counted as single particles, the expression 
of Smoluchowski for the total number of ag gregates is recovered. This result was repeatedly 
validated by ESS of non-deformable drops] 1031 104 1. In fact, an alternative expression for the 
total number of aggregates was deduced 1 103ll . This expression is valid for systems of distinct 
polydispersity with 0.01 < (p < 0.30, 10~ 5 M < C sur f < 8 xlO 5 M (where C sur f is the surfactant 
concentration). In the absence of a repulsive barrier, limiting values for the mixed flocculation- 
coalescence rate (kfc) of non-deformable drops were also calculated II 1 0311 - The value of kfc 
depends on the volume fraction of oil and the Hamaker constant (Ah) of the van der Waals po- 
tential between the drops. At infinite dilution kfc is equal to 6.4 xl0~ 18 rrr's^ 1 for Ah = 1-24 
xl0~ 21 J, and 9.5 XlO -18 wzV 1 fovA H = 1.2 xl0~ 19 J. 

The aggre gatio n d ynamics o f non- deformable [To! 0, CM CM H and de- 



formable 1561, 11071 |6j U081 106] drops subject to coalescence has been extensively studied by 



our group using simulations. Since the scale of drops and surfactants differs in several orders of 
magnitude, the explicit movement of both entities at the same time is not yet computationally 
feasible. Hence, the major problem of ESS is to consider surfactant effects without moving the 
surfactant molecules explicitly. The second major difficulty is to account for the influence of a fi- 
nite surfactant population in all destabilization processes simultaneously. The surfactant surface 
excess depends on time due to the continuous decrease of the available interfacial area caused by 
the coalescence of the drops, and the molecular diffusion of the amphiphile lfl09ll . The adsorp- 
tion of surfactant determines the interfacial properties of the drops such as tension, charge, fluid 



mobility, etc. These properties change the deformability of the drops, as well as their interac- 
tion potentials, hydrodynamic tensors, and the coalescence mechanism (film drainage, capillary 
waves, hole formation) (^(J. Hence, all destabilization processes depend on the surfactant distri- 
bution throughout the system, and this distribution is continuously changing. 

Former ESS calculations were exclusively devoted to non-deformable drops and ionic sur- 
factants. Although ESS confirmed most of the predictions of the DLVO theory, some important 
differences appeared. According to the simulations of nondeformable drops, the stability thresh- 
old occurs as soon as the maximum of the repulsive potential starts to form IU06II : j^- = and 

= 0. The behavior of the systems depends strongly on the surfactant concentration and the 

total interfacial area of the emulsion. If the surface excess of the surfactant is approximately 
constant, as it occurs for very low or very high surfactant concentrations, DLCA (Diffusion 
Limited Cluster Aggregation) and RLCA (Reaction Limited Cluster Aggregation) regimes are 
observed, respectively II96I1 . However, if the redistribution of surfactant molecules after a coa- 
lescence event is considered, the surface potential of the drops evolves dynamically, promoting 
a progressive stabilization of the emulsion with respect to flocculation and coalescence jl06ll . 
In these cases, the fractal coefficient of the aggregates changes with time, and the equation of 
Smoluchowski has to be modified to account for the total number of aggregates. Some recent ex- 
periment al evidence on the existe nce of surfactant redistribution upon coalescence can be found 
inRefs. iHTol [ml [TT1 [TT1 [TT1 . 

Using silicon-in-water nanoemulsions of 276, 374 and 576 nm and SDS concentrations of 
3, 5, and 7 mM, Barchini and Saville ll82ll showed that the zeta potential predicted by mobility 
and conductivity measurements cannot be reconciled. The system of 576 nm with 5 mM even 
shows a negative slope of the conductivity as a function of the volume fraction of oil (<p < 0.05). 
These inconsistencies were ascribed to the inadequacy of the standard electrokinetic theory to 
describe the high surface potentials of the droplets and the properties of the inner part of the 
diffuse double layer. However, recent evidence indicates that the interfacial area occupied by the 
surfactant in nanoemulsions could be one order of magnitude higher than the one expected from 
the equilibrium adsorption of the surfactant to macroscopic interfaces 111 1511 . 

Poulin et al. Ill 1611 noticed that the behavior of silicon-in-water nanoemulsions (R = 190 
nm, [NaCl] = 0.35, [SDS] = 2mM) depends on the temperature and on the volume fraction of 
oil. Following the changes of the scattering peak, the aggregation mechanism was identified as 
DLCA for deep quenches and a volume fraction of oil between 0.01 < (f> < 0.1. For <p > 0.1 and 
shallow quenches, different mechanisms closer to spinodal decomposition were observed. It was 
also noticed that the droplets become highly adhesive at high salinity exhibiting large contact 
angles. According to Wang et al. II 1711 . flocculation is the main mechanism of destabilization in 
water/decane emulsions stabilized with a mixture of pentaoxyethylene lauryl ether (C12E5) and 
di-dodecyldimethylammonium bromide. 

As any o/w em ulsion, nanoemulsions are subject to creaming and Ostwald ripening [118, 



1 19l 12d 121 , 122 1. However, nano-metric particles should not be appreciably subject to cream 



ing unless they aggregate. In the absence of salt, all systems investigated showed an appreciable 
amount of creaming after two days. Creaming might be prevented mixing dodecane with a high- 
density oil, to make the mixture neutrally buoyant. Unfortunately, most high density liquids are 
halogenated alkanes, whose solubility in water is higher than the one of dodecane. This favors the 
occurrence of Ostwald ripening. The theoretical ripening rate of surfactant-free dodecane/water 
emulsions is equal to 1.3 xl0~ 25 m 3 s~'. However, this rate changes with surfac tant a ddition. 
Ostwald ripening might be prevented mixing dodecane with an insoluble oil II 12311 12411 . A sys- 
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tematic study on the influence of creaming and ripening on the value of the flocculation rate is 
currently carried on in our lab. In this article we concentrate on the evaluation of aggregation in 
a typical nanoemulsion. 

It is evident from above that emulsions are far more complex than suspensions. As a result, 
the quantification of the processes leading to destabilization is more difficult, even if the rele- 
vant variables involved are clearly identified. The predi ction of the t emporal evolution of these 
systems requires elaborate simulation techniques ill 06 . 72 . 56l 107 1. and the quantification of 
each destabilization process is restricted to those conditions in which the rate of each mechanism 
is substantially different from the rest. Moreover, it is uncertain whether these rates are rele- 
vant in typical systems where all processes occur simultaneously, regulated by a time-dependent 
surfactant distribution. 

In this work, the initial 30-seconds variation of the turbidity of an ionic dodecane-in-water 
nanoemulsion is used to evaluate its flocculation rate at different ionic strengths. This time span 
is considerably larger than the one required for the formation of d oublets. The na noemulsions 
studied were prepared by the phase inversion composition method 1 125 , 126l 127 , 128 1. which 
basically consists in the sudden dilution of a pre-equilibrated system in the balance zone. The 
experimental determination of W as a function of the initial slope of the turbidity (r) vs. time (f), 
is complemented with four theoretical models of the turbidity which depend on the average drop 
size and the flocculation rate. These models are approximate formalisms which allow quantifying 
the stability of nanoemulsions with respect to flocculation and coalescence. Adjustment of the 
theoretical expressions to the experimental data allows the evaluation of either a flocculation rate 
or a mixed flocculation-coalescence rate. This methodology provides a complementary method 
to evaluate W, which proves to be in agreement with direct experimental determinations. 

In the following section, the necessary theoretical background of the work is presented. In 
the third section the experimental procedure is outlined. The fourth section contains the results 
of the measurements and the theoretical fittings. In the fifth section, the results are discussed. 
The paper finishes with some concluding remarks. 



2. Theoretical background 

According to the scheme proposed by Smoluchowski, the process of aggregation begins with 
the formation of doublets: 

dvii 2 

— = -k u n x - k n n\n 2 - ... (5) 

dn 2 fen 2 , IC \ 

— = — Tij - k\ 2 n\n 2 - ... (6) 

where n\ and n 2 are the number concentration of single and double particles, and k\ \ and k\ 2 are 
rate constants. 



In general, the turbidity of a colloidal dispersion can be defined as 1 26, 28]: 

oo 



T ■ 

k=\ 



where is the number of aggregates of k particles per unit of volume, and cr> is their total light 
scattering cross section. Differentiation of Eq. Qgives: 

dr dn\ ^ dn 2 ^ 

dt dt dt 
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For the limiting case of t — > 0, n% approaches zero and the second and further terms on the right 
side of Eqs. [5]|6l and[7]can be disregarded. Combining Eqs. [5] [6] [7] and[8] 

(f)-°(^Hr>-4..< 

Here, o~\ and <T2 are the optical cross-sections of a single spherical particle and a doublet, t = 
(l/L)ln(/o/^), Abs = log(/o/7) is the absorbance of the dispersion, / is the intensity of light 
emerging from a cell of width L (10~ 2 m), and Io is the intensity of the incident light. 

The above cross sections can be evaluated using a suitable theory of light scattering. In this 
case, k\\ can be deduced from turbidity measurements. 

Rayleigh, Gans, and Debye [ 35] formulated an approximate theory of light scattering, re- 
stricted to particles with small relative refractive index, and small size compared to the wave- 
length of the incident light. This theory is valid whenever: 

C RGD = (AnalA)(m~Y)<&\. (10) 

Here, A is the wavelength of light in the medium, and m is the relative refractive index between 
the particle and the medium. 

Kerker 13511 estimated the regions of Crcd vs. a for which the error in the extinction cross 
section of a sphere is lower than 10%. For latex particles (m =1.20), a must be lower than 2.0. In 
the case of dodecane (m =1.07, A =400nm, a — R =72.5 nm) a =1.51, and Crcd =0.22. Hence, 
the errors produced in the calculation of the optical cross section of the drops using RGD are 
expected to be smaller. 

W can be obtained from the initial slope of the curve t v s. / at low and high ionic strength 
(slow and fast aggregation, respectively) 1 129 , 130i 131 , 132 . 27. 251: 



w _ w (dT/dt) QJasl _ (dAbs/dt\ fast 

exp {dTldt\ slow (dAbs/dt\ slow ' 

Because in the case of solid particles the optical factors of Eq. |9]are canceled yielding: 

k fast 

W = W *n = 7^- d2) 

It is also possible to fit the experimental data on the total number of aggregates to Eq. [2] to 
obtain the value kf. In this case, the stability ratio can be written in terms of the values of the 
flocculation rate in the absence (k^ ast ) or presence (k s J° w ) of a repulsive force: 

W=-t- . (13) 

KSlOW 

K f 

Thus, if the flocculation rates of aggregates of different size are similar, the average effect of the 

uslow _ ufa 

V ~ r / 



repulsive barrier can be incorporated in Eq. ^substituting kf by k s l ow = l/, as ' /W. 



2.1. Theoretical models for the estimation of the turbidity curve 
2.1.1. Aggregation of emulsion drops (model Ml) 



Puertas et al. ll 1 3 3 . 134 1 developed a general method for evaluating the average cross-section 
of an aggregate of size k. The total cross-section of each aggregate results from four contribu- 
tions which take into account the possible structures that might be generated when each pair of 
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particles in an aggregate is separated by zero, one or two particles between them. Within the 
R GD ap proximation, the cross-section of an aggregate of k spheres of radius a, can be expressed 
as lfl35lf35ll: 

4 r" 

cr k = -na 2 a\m - l) 2 P k (ff)(l + cos 2 ■&) sia(ff)d§. (14) 



with: u = 2a-sin(#/2), and 



Here § is the scattering angle, and P k (&) is defined as: 

P k (ff) = P 1 (&)(k + A k ), (15) 

where Pi is given by: 

/ sin m - m cos « \ 2 

pm = (3 — - 3 — ) , (i6) 

Here, s, 7 is the center-to-center distance between spheres i and j, and h - u/a = {An/ A)sin(fl/2). 

For the scattering cross-section of aggregates with k > 4, Puertas et al. considered in the 
above summation only the terms corresponding to pairs of particles separated by two other par- 
ticles at most, because they decrease with the distance between them (note that pairs of particles 
close to each other but not connected by two particles or less are thus not accounted, which 
restricts the validity of the calculation to small aggregates). Following this procedure, the aver- 
age value of A k can be calculated if it is considered that all structures of k particles are equally 
probable. In this case: 

< k> - x 

i,r, i*j 1 

2(^1)^,2(^2)1 r^sin ( 4 M co s(y / 2) ) 
2ha 2n Jo 2u cos(y/2) 

16n 2 J_ 2jr/3 J-2„/3 hl(y u y 2 ) 

with , J2) equal to: 

2 (n + 72) 



'(71.72) = 2a 



1+2 cos(yi) + 2 cos(y2) + 4 cos 



(19) 



Angles y, j\ and 72 are shown in Fig. Q] 

Eqs. [18] and [19] can also be calculated using an alternative set of equations which is formally 
equivalent, but more computationally efficient. These equations are given in Appendix A for 
further reference. Using the expressions of Smoluchowski for the kinetics of aggregation, the 
turbidity can be written as: 

00 

tmi =2^n k {kf,t)(r kia , (20) 

k=\ 

where cr^a = o~k is the optical cross section of the aggregates, and: 
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no(k f n Q t) k 1 

nt = j—r- (21) 

(l+k f n t) k+1 

Puertas and de las Nieves IU33I1 used the kinetic constant, and the time to as fitting parameters. 
Notice that to is the time at which the process of aggregation begins. Theoretical calculations 
were compared with experimental evaluations of the aggregation of latex suspensions. A good 
agreement between theory and experiment was found. The theoretical approach reproduced 
the experimental data of t corresponding to 50 seconds of aggregation. The kinetic constants 
obtained agreed with the ones calculated from the initial slope of t vs. t, and the literature 
(fen = 2k f), indicating that the kinetics of Smoluchowski is a good approximation to the real 
aggregation process. 

Thus, the model Ml refers to the direct application of t his m ethodology to the case of emul- 



sions. Despite the fact that the model developed in Refs. 1 133, 134] describes the aggregation 
behavior of a monodisperse system, we aim to apply it to a more real system, where polydisper- 
sity is a natural feature. It is clear that any distribution of particle size has a certain degree of 
polydispersity. However, the theory of Smoluchowski addresses the case in which all particles 
have the same radius at the beginning of the process. The fact that this formalism is commonly 
validated using suspensions (despite the inherent polydispersity of the systems) is an implicit 
recognition that the optical techniques accept a certain degree of tolerance in the particle size. 

2.1.2. Complete coalescence of the drops after contact (model M2) 

Model M2 assumes that the drops do not aggregate, but instead coalesce as soon as they 
collide with one another. 

Drops coalesce if either the attractive potential or the thermal interaction with the solvent 
provides enough energy to overcome their repulsive interaction. In this case, a new particle with 
radius at is created at the center of mass of the coalescing drops. Using the conservation of 
volume, a 3 + a 3 = a 3 , hence, 

a k = fka. (22) 
Since, the form factor of a single sphere (Pk, s ) of a radius a k is defined as: 

sin u k - u k cos u k 
rk, s (y) = 5 



(23) 



The cross-section of a sphere of radius a k (cr^.,) is given by: 



4 c 

o- k , s = -na\a\(m - l) 2 P k J$)(l + cos 2 ff) sin §d§, (24) 
9 Jo 

where a k = 2na k /A, and u k - 2a k sin(#/2). 

If the drops do not deform, but coalesce as soon as they flocculate in the primary minimum, 
the turbidity can be calculated as: 

oo 

tmi = 2 n k (k f , f)cr kyS . (25) 

k=\ 

This model is meant to reproduce the cases in which there is neither a barrier for flocculation 
or coalescence. Previous emulsion stability simulation demonstrated that in this case, the dy- 
nam ics o f Smoluchowski applies, since the number of aggregates change as the number of drops 



II961. 110611 . Notice however, that the number of particles decreases as a product of coalescence. 
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2.1.3. Partial coalescence of flocculating drops ( model M3 ) 

An initially monodisperse emulsion rapidly turns into a polydisperse system as a conse- 
quence of the fiocculation and the coalescence of the drops. The turbidity of this system results 
from several contributions including primary (small) drops, aggregates of primary drops, bigger 
spherical drops resulting from the coalescence of small drops and "mixed" aggregates resulting 
from either the fiocculation of big and small drops, or the partial coalescence of aggregates of 
small drops. If the number of "mixed" aggregates is negligible (see Appendix B): 

Kmx Kux 

t M 3 = nicri + 2^ n Ka cr Ka + ^ n Ks o~ Ks , (26) 

k=2 jt=2 

where k max represents the maximum number of particles in an aggregate of size k. Subscripts a 
and s stand for aggregates and spherical drops, respectively. 

If the large drops are originated during the process of fiocculation of smaller drops, their 
number should be proportional to the values of n k proposed by Smoluchowski. The number of 
collisions predicted by Smoluchowski will be the same, but only a fraction x k jfl of those collisions 
will produce aggregates of size k: 

rik. a = n k x kA . (27) 

If it is assumed that the fractions of aggregates (x kya ) and drops (x ks ) are the same independent 
of the value of k, then: 

rik, a = n k x Ka = gn k , (28) 
n k , s = n k x Ks = (1 - g)n k . (29) 

Therefore, 

kmx kmx 

t m3 = mo-y +g^n k cr Ka + (1 - g) ^ n k cr kiS . (30) 

k=2 k=2 

In the theory of Smoluchowski, an aggregate of k particles contains k primary particles. 
Therefore, its volume is equal to k times the volume of a single primary particle. Hence, if 
the process of coalescence does not occur during the fiocculation of drops but afterward, it will 
still affect the number of aggregates n k in the same way. The fraction x ktS results from both 
the instantaneous coalescences of the drops at contact or their "delayed" coalescences within an 
aggregate previously formed. 

When coalescence does not occur and only the mechanism of aggregation operates, the frac- 
tion of floes produced by aggregation is equal to one {x kt „ = g = 1). In this case, the third term 
on the right hand side of Eq. [30]disappears, and Eq. |20]is recovered. 

Notice also that the summations of Eq. [30]can be regrouped, giving rise to an average cross 
section for floes and drops: 

k„ax 

t m3 = mo-\ +^n k [g o- Ka + (1 - g) o- ktS \. (31) 

k=2 

2.1.4. Modification of Ml to account for a finite coalescence time (model M4) 

Model M2 supposes that the coalescence time is negligible in comparison to the time required 
by the incident laser light to sample the structures of the system. If that is not the case, it might 
be possible to "observe" structures in which the flocculated drops are deformed but have not 
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totally coalesced. This consideration is specially important if the deformation is such that (an 
approximately planar) film is formed between the drops. Small drops are not expected to deform 
significantly prior to coalescence. However, they have to form some sort of "bridge" between 
their internal phases before they merge into one another. In this case, the form factor of the 
transient structure would be intermediate between that of two flocculated particles and a bigger 
sphere. In order to test this possibility, we introduced an average cross-section for the transient 
structure equal to: 

Ctransl - ~ " , (.J A) 

(1 - p)cr Ka + p (T k<s 

where p is a real number and < p < 1. Notice that the use of equations of the type: 

O- trans! = P 0~kM + (I - p) <T Ks , (33) 

leads to the same expression of turbidity resulting from model M3. 
Use of Eq. [32] allows to express t in the following form: 

r M4 = ni(T\ + V n k cr mmsl . (34) 

k=2 

This equation is equal to Eq. [3T]except for the expression of the cross section used in this case: 
o~trans\ - Thus, when p — only the spherical drops will contribute. In this case, the optical cross 
section cr rra „ sl (Eq. [32b becomes equal to cr^ s and Eq. [25]is recovered. This represents the case 
where only coalescence operates. Conversely, if only flocculation occurs p = 1, cr trans \ = 0> a 
and Eq. [20]is recovered. 



3. Experimental 

3.1. Materials 

N-dodecane (Merck, 98% purity) was purified by passing it twice through a column with 
alumina. Sodium dodecyl sulfate, SDS (Merck) was recrystallized from ethanol two times. NaCl 
(Merck, 99.5% purity) and isopentanol (Scharlau Chemie S. A., 99% purity) were used as re- 
ceived. The water used in all experiments was distilled and deionized (conductivity < 1/iScm _1 
at 25 °C) using a Simplicity purificator (Millipore, USA). 

3.2. Phase diagram [NaCl] vs. [SDS] 

In order to select a starting system for nanoemulsion preparation, unidimensional scans of 
salt concentration (wt% NaCl) were made at several surfactant concentrations (wt% SDS). The 
weight-fraction of water was kept constant at f K = 0.20 {(p = 0.84). The isopentanol concentra- 
tion was also maintained at 6.5 wt.%. This alcohol was added in order to facilitate the transfer 
of the surfactant to the oil phase. Each system was prepared in vials, weighing each compo- 
nent to the final composition. The substances were added in the following order: water, SDS, 
isopentanol, dodecane and NaCl. The vials were subsequently gently shaken with a Varimix 
(Barnstead International, USA) for 2 hours in order to facilitate the contact between phases. 
Finally, the systems were placed in a thermostated bath at 25 °C until equilibrium was reached, 
i.e., until complete phase separation was observed. The number of phases was evaluated visually. 
Anisotropic liquid crystalline phases were identified using polarized light. 
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3.3. Nanoemulsion preparation 

Water/NaCl/isopentan ol/SDS/dod ecane nanoemulsions were prepared using the phase inver- 
sion composition method [ 128. 1127I1 . AW, + LC + O system was selected as a starting system. 
This system was suddenly diluted at constant stirring until a final weight-fraction of water of 
f w = 0.95 (<p = 6.5 x 10~ 2 ). This procedure allowed to obtain a concentrated mother nanoemul- 
sion. The final nanoemulsions were obtained diluting the mother nanoemulsion with a suitable 
aqueous solution of SDS. During the dilution, the surfactant concentration and the number of 
drops per unit volume were adjusted to 8 x 10~ 3 M and 2.0 x 10 11 drops/mL, respectively. The 
number density of the drops was estimated considering the mean radius of droplets obtained from 
photon correlation spectroscopy (PCS) measurements. 

3.4. Nanoemulsion characterization 

The ^"-potential of the drops was evaluated using a Zetasizer-Z (Malvern, UK). This equip- 
ment employs the Smoluchowski approximation IU36I1 which relates the mobility to ^-potential. 
The average droplet size and the droplet size distributions (DSD) of the nanoemulsions were 
evaluated by PCS using a Malvern Zetasizer 2000 (Malvern, UK). 

The form of the DSD at a selected set of times was determined using a BI-200SM goniometer 
(Brookhaven Instruments, USA). 



3.5. Experimental determination of stability ratio (W) 

In order to determine the stability ratio of a nanoemulsion, the change of the absorbance as a 
function of time at different electrolyte concentrations was measured. The optimal experimental 
conditions were obtained after several preliminary studies. These included: a) the adjustment of 
the velocity of addition of the salt solution in order to ensure the complete homogenization of 
the sample, and b) the determination of the optimal concentration of drops (from dAb s/dt vs. 
[Drops/mL] curves at different wavelengths) in order to avoid multiple scattering lll37ll . 

A Genesys 5 spectrophotometer (Thermo Scientific, USA) was used, at a wavelength of A — 
400 nm. At this lambda, the emulsion components do not absorb significantly, and the fulfill- 
ment of RGD conditions is guaranteed. The methodology employe d was very similar t o the one 
previously used for studying the aggregation of latex suspensions Il34lll38lll32lll3lll : 2.4 mL 
of the nanoemulsion was placed into the cell of the spectrophotometer and its optical absorbance 
was measured. Following, 0.6 mL of the respective NaCl solution was quickly added in order to 
reach the desired electrolyte concentration. The salt concentration was changed between 380 and 
600 mM. In all cases, the concentration of surfactant in the cell was kept constant at 8 x 10~ 3 M. 
Higher surfactant concentrations promote the formation of hydrated crystals of SDS above 600 
mM NaCl. 

The optical absorbance as a function of time was recorded continuously using a computer 
for a period of 30 seconds. The measurements were repeated at least three times at each salt 
concentration. Since the aggregation process is very fast, the initial value of the absorbance 
after the addition of salt solution is unknown. An approximate value of the initial absorbance 
(Abs([NaCl = 0])) was obtained, diluting 2.4 mL of the nanoemulsion with 0.6 mL of pure water. 
This was done previously to each experiment. From the value of this absorbance, the initial time 
of aggregation (to texp ) was estimated in the actual experiment. 
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3.6. Evaluation of W exp 

As a first approximation to W, th e initial slope o f the curves of Abs vs. t was considered to be 
proportional to the coagulation rate 1 3 3l 1 1 3 ill . Notice that neglecting the difference in the 



radius of different preparations of the nanoemulsion, the cross-sections of the aggregates can be 
simplified in Eq. |9]giving in Eqs. QT|and[T2] The value of W exp (Eq. [TO was approximated by the 
quotient between the average initial slope obtained under fast aggregation conditions (between 
420 and 600 mM of NaCl), and the one obtained at diluter electrolyte concentrations (380-410 
mM NaCl). 

3. 7. Evaluation ofWk n 

Using Eq. [9] and the initial slopes of Abs vs t at each salt concentration, the rate for doublet 
constant formation (kn) was calculated. The cross-sections, en and cr 2 were estimated using Eq. 
[14] for k = 1 and k — 2, respectively. These equations include the average radius (R exp ) of the 
emulsions. Since one mother nanoemulsion was prepared for each salt concentration, the values 
of en and a"2 are slightly different for each preparation. Hence, the quotient Wt n = k^"'" /kf° w 
(Eq. fT2l is not necessarily equal to W exp . 

(, \ slow 

-cr 2 -o-\) 

Wexp — -r- = W ku (35) 

exp , \fast "11 v ' 

3.8. Fitting of the turbidity curve 

The models Ml, M2, M3 and M4 were programmed using a symbolic manipulation program 
(Mathematica 8.0.1.0). This program uses least-squares fittings in order to adjust the theoretical 
model of the turbidity to the experimental data. In all cases, the absorbance data was converted 
to turbidity using the relation: t = 230Abs. The values of kf, to, g, and p are obtained from the 
fittings. The code of model Ml was tested using the original data from Ref. lll33ll . reproducing 
the same kinetic constant previously obtained. 

The average radius (R /,-,) of the drops is a parameter of the calculation, which in some cases 
was slightly varied in order to: a) maximize the quality of the fitting, and b) obtain a theoretical 
value of to close to the experimental one (to^ xp ). 

Once the flocculation rates were obtained from the theoretical models, the stability ratios 
were calculated for each salt concentration. 



3.9. Polydispersity 

In order to appraise the effect of the initial polydispersity of the emulsions on the quality of 
the fittings of kf, we recalculated the flocculation rate using the normalized frequency of the drop 
size distribution for three particle radii: R m j„ = Rfi, - a s d o~ s d, and R max = Rfu + a s dO~sd, where 
cr s d is the standard deviation of the distribution and a s d = 1,2. Using these radii, an approximate 
value of the turbidity can be calculated as: 

T aprox = T(R m i„)f Rwm + T(R fil )f Rfjl + T(R max )f Rwax , (36) 

where fit stands for the frequency of the DSD corresponding to radius R (/«„„„ + fii fu + fg^ = 1). 
The effect of and Rmax was observed, comparing the behavior of T aprox with the one of t(R /,,) 
for models Ml and M3 at [NaCl]=390 mM and 475 mM. The difference between R exp and R fit 
for 390 and 475 mM NaCl is 14 and 17 nm, respectively. 
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4. Results 



Fig. [2] shows the phase behavior of the systems as a function of surfactant and salt con- 
centrations. The dash lines indicate the approximate limits of different phase behaviors. The 
three-phase region is located in the center of the diagram. This region is divided into two zones: 
a) a subregion containing water (W,), a bicontinuous microemulsion (D) and an oil phase (O), 
which is formed at a low surfactant concentration (4 wt.% SDS); and b) a second subregion 
composed of water (W,), liquid crystal (LC) and an oil phase (O), formed at high surfactant con- 
centrations (> 6 wt.% SDS) 1139;]. Below the three-phase region, a two-phase region consisting 
in a direct micellar solution (W m ) and an excess oil phase (O) is observed. Here, the surfactant is 
preferentially dissolved into the aqueous phase and the formation of o/w emulsions is favored. 
Above three-phase zone, a two-phase region exists, consisting of excess of water (W t ) in equilib- 
rium with a reverse micellar solution (O m ). When the surfactant concentration is higher than 6 
wt.%, a one-phase region appears, O m . 

As Fig. |2] indicates, it is possible to promote a phase transition changing the salt concen- 
tration. During a transition involving three phases, the system shows low values of interfacial 
tension, which favo rs the formation of drops of nanometer size. According to previous experi- 
mental works jl40l 141 , 142 , 143 , 126 1. the presence of liquid crystals favors the formation of 
stable nanoemulsions. Hence, Water/NaCl/isopentanol/SDS/dodecane nanoemulsions were pre- 
pared starting from a W t + LC + O system. 

The size distrib ution of a typical nanoemulsion prepared by the phase inversion composition 
method jl28lll27ll is shown in Fig. [3] It has a log-normal distribution with an average radius 
around 73 nm. The physical characteristics of the dispersion reported in Table [2] correspond 
to the average of ten nanoemulsions prepared independently. The composition of the starting 
system, the concentrated mother emulsion and the final diluted nanoemulsion are shown in Table 
1. 

Previous to the turbidity measurements, the evolution of the average droplet radius (R exp ) of 
several nanoemulsions as a function of time was studied (see Fig. @J. In the case of emulsions, 
it is uncertain whether the fluctuations observed in Fig. ^correspond to the statistical error of 
the measurements or they reflect a true variation of the average radius. In any event, it is clear 
that the average radius does not change considerably within the first 30 minutes. This time it is 
substantially longer than the one required for the evaluation of the turbidity at each ionic strength. 
Nevertheless, one mother nanoemulsion was synthesized for each salt concentration. 



4.1. Variation of turbidity as a function of time 

Fig. [5] shows the typical curve of Abs vs. t of an aggregating colloidal dispersion. Prior 
to addition of the salt solution, the absorbance is constant. As soon as the 0.6 mL of the NaCl 
solution is added, the absorbance decreases due to the dilution of the emulsion. Then, it increases 
indicating the aggregation of the drops. As can be observed in Fig. [3] there is a non-monotonous 
variation between 4 and 6 seconds. During this time, two events occur: a) the homogenization 
of the sample due to the addition of salt solution, and b) the undesired input of external light that 
enters the measuring cell when the compartment of the instrument is opened in order to add the 
salt solution. The process of flocculation begins at time to <exp which occurs during this transient 
period. This value is used as a reference for the theoretical value, toju, obtained from the fitting. 
The comparison of these two values is an additional test of consistency of the theoretical models. 

The relevant variation of Abs vs. t begins as soon as the monotonous curve is observed. At 
this point, the initial slope of the curve can be calculated. 

18 



Fig. [6]shows the initial slope of Abs vs. t as a function of the salt concentration. Between 380 
and 420 mM NaCl, the emulsion is slowly destabilized. In suspensions, this regime is associated 
to Reaction Limited Cluster Aggregation (RLCA). Above 430 mM NaCl, the emulsion destabi- 
lizes very fast, showing values of dAbs/dt which are similar within the error bars (Appendix C). 
This regime corresponds -in the case of suspensions- to Diffusion Limited Cluster Aggregation 
(DLCA). 

Fig. [7] shows the results of two measurements of the initial droplet size distribution of the 
same nanoemulsion using the algorithms of Contin jl44ll and Cumulants analyses jl45ll imple- 
mented in the Brookhaven BI-200SM Goniometer. Fig. [8] illustrates the evolution of the DSD 
for an "RLCA" system ([NaCl] = 380 mM). The average radius and the width of the distribution 
increase monotonically as a function of time due to the aggregation of the drops. The width of 
DSD enlarges considerably after 600 s. At high ionic strength (Fig. |9l [NaCl] = 600 mM), the 
polydispersity of the system increases much more pronouncedly due to the coalescence of the 
drops (notice the scale of the abscissa). The distribution turns bimodal after 60 s, with peaks 
around 400 and 2000 nm according to the Contin analyses. Between these two extreme cases, a 
certain degree of stabilization is exhibited (Fig. |6}. Notice the difference between the DSD cor- 
responding to 395 and 400 mM (Fig. ITOb . This i s the coalescence threshold suggested by models 
M3 and M4. According to the simulations lll04ll . the progressive decrease of interfacial area due 
to the coalescence of the drops and the dynamic adsorption equilibrium of the surfactant to the 
available interfaces generate a progressive increase of the repulsive forces between the drops as 
a function of time. If the drops remain non-deformable, the total surfactant concentration and 
the ionic strength of the solution will determine if complete phase separation occurs or a lower 
value of the total interfacial area of the emulsion can be stabilized. 

4.2. Fitting of the turbidity curve with theoretical models 

Table[3]shows the kinetic constants obtained for doublet formation {k\\). Figs. ITT1 IT2l and [T3l 
show a comparison of the theoretical models for t with the experimental measurements for the 
cases of slow (Figs. [TT1and[T2l and fast aggregation (Fig. [T3l . 

The data in Fig. QT]correspond to [NaCl] = 390 mM. The theoretical models were calculated 
using k max = 100 (Appendix D). While model Ml reproduces the data quite well, the model of 
spheres (M2) roughly follows the experimental trend. Fig. [T2]compares the same experimental 
data to the predictions of models M3 and M4. Both models fit the data satisfactorily. However, 
they showed a slight deviation at the initial times. Notice that the average aggregation time 
tave = h/2 — 1 /nokf is of the order of 9.3 s. Hence, the models are able to reproduce at least 3 

In the case of fast aggregation (Fig. [T3l . [NaCl] = 600 mM, models Ml and M2 fail to repro- 
duce the experimental data completely, while the models M3 and M4 still hold. It is impossible 
to distinguish between M3 and M4 in Fig. [T3]since both fittings are of similar quality. 

For completeness, Fig. [14] shows the turbidity curves calculated with model Ml for k max = 
2, 3, 4 and 100, even for this 30-seconds evolution. As expected a large value of k max is necessary 
in order to reproduce the experimental curve completely, due to the high rate of the reaction. 

We move now to comparing the stability ratios. W exp and Wk u were obtained using Eq. QT| 
and Eq. Q~2] respectively. A plot of W exp vs. [NaCl] clearly shows the regions corresponding to 
the slow and fast aggregation (see Fig. IT3T l. Using the kinetic constants of Table [3] the values 
of Wjt n were calculated and compared with the experimental ones (W exp ). There is a very good 
agreement between both determinations. 



19 



Tables [4] and [6] show the kinetic rates of fiocculation obtained from the fitting of models 
Ml, M3 and M4. In Table|H the kinetic constants between 380 and 395 mM NaCl result from 
an appropriate fitting of the data by the theoretical model (Ml). Above 400 mM NaCl, the 
model does not adjust to the experimental trend. Consequently, the values of the error bars of the 
kinetic constants are higher than the average value. Note that the particle radius, that has been 
chosen as a fitting parameter (Tables [5] and |6) fluctuates around the value measured with PCS 
(Fig. |4j, showing that it can be taken as a robust measure, despite the inherent polydispersity of 
the system. Interestingly, this value of the particle radius is obtained from a fitting of the whole 
turbidity curve, i.e. where large aggregates are also included, and not only single particle data. 

The values of W were determined from the data of Tables [5]and|6] Fig. [l6]shows a compar- 
ison between the theoretical and experimental values of the stability ratio. The fast fiocculation 
rates required to calculate Wk n , Wmt, and Wma were approximated by the average of the values 
of k f obtained between 420 and 600 mM NaCl (kf™ =3.2 xl0~ 18 mV 1 , k^" =1.3 xl0~ 18 
m 3 * -1 , and k{^f = 1.2 xl0~ 18 m 3 s _1 ). The agreement is satisfactory considering the complexity 
of the emulsion system. In the case of emulsions, the turbidity is associated with the change in 
concentration of the droplets, which is due to the phenomena of fiocculation, coalescence and 
Ostwald ripening. While in the case of solid suspensions, the initial slope of turbidity curve is 
only due to aggregation, and more specifically, to the formation of doublets 12711 . 

The values of the critical coagulation concentration (CCC) of the nanoemulsion predicted 
from the experimental data and the theoretical fittings (Fig. [T6l > are shown in Table [7] 

4.3. Effect of the initial polydispersity on the evaluation ofkf 

Figs. fT71 and [TBI illustrate the effect of the initial polydispersity on the calculation of kj. As 
the polydispersity increases, the theoretical value of the turbidity deviates from the experimental 
points. This happens for both models Ml and M3 at both salinities tested (390 and 475 mM 
NaCl). However, the deviations are larger in the case of model Ml at 475 mM since it does not 
adjust the experimental data. Since model M3 converges to model Ml when only aggregation 
occurs, the deviations of these models are very similar for [NaCl]= 390 mM. Similar results 
are obtained using model M4 instead of M3. In both cases (M3 and M4) an increase in the 
polidispersity of the system provokes a displacement of the turbidity curve, without changing its 
curvature. This implies that a smaller radius should be used to reproduce the experimental data, 
but leaves our main conclusion unchanged, namely, the interplay of coagulation and coalescence 
have to be considered to describe correctly the destabilization of nanoemulsions, such as is done 
in our models M3 and M4. 

5. Discussion 

As shown in Fig. [6] a dodecane-in-water nanoemulsions stabilized with SDS shows the 
typical behavior of charged latex suspensions with respect to the salt concentration. Below the 
critical coagulation concentration (CCC =411 mM NaCl) the system is kinetically stabilized 
with respect to fiocculation. The stability decreases progressively with the increase of the ionic 
strength. Above the CCC, the fiocculation of the particles occurs very fast, and the aggregation 
rate reaches a plateau, similar to the one exhibited by solid particles under DLCA conditions. 

It is important to remark that each point in Fig. [6] corresponds to the average at least three 
different measurements. During the first 30 minutes, the emulsion is fairly stable: its average 
radius changes less than 10%. After this time, the average radius changes considerably even 
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in the absence of an appreciable concentration of salt. This instability probably occurs as a 
consequence of Ostwald ripening IT72lll40ll . which is favored due to the high polydispersity of 
the emulsions and the appreciable solubility of the dodecane in water. 

In the case of suspensions, the addition of salt screens the surface charges of the particles 
promoting their flocculation. In the case of emulsions, the charge of the surfactant adsorbed to 
the oil/water interface is also screened, but further surfactant adsorption is also promoted. Hence, 
the stability of the emulsion results from a balance between these two phenomena. Once the salt 
concentration is high enough to produce a maximum surfactant adsorption, subsequent addition 
of salt will only generate the progressive screening of the surface charges as predicted by the 
theory ofDLVO. 

At [NaCl] < CCC (see Table [7]), the original drops (primary particles) are expected to ag- 
gregate following the dynamics proposed by Smoluchowski Ir96lll06ll . Thus, Eqs. l2l[T3landl2TI 
should be observed. The initial slope of the curve of Abs vs. t at each salt concentration can be 
used along with Eq. [9] to obtain the rate of doublet formation {k\\). Table [3] shows a minimum 
value of 2.1 x 1CT 19 m 3 s~ l . The highest value obtained corresponds t o 3.8 X 10" 18 m 3 s _1 . This 
value is very similar to the one reported by Puertas and de las Nieves ll 1 3 3h for the case of latex 
particles (MP3, R = 92.5 nm). As shown in Fig. Q3] the stability ratios W exp and Wit,, are very 
similar. 

Model Ml represents the effect of the aggregation of the drops on the turbidity. This model 
fits the data for 380 < [NaCl] < 395 mM NaCl within the slow aggregation regime. This is 
illustrated in Fig. QT|for the case of [NaCl] = 390 mM. As shown in Fig. [T4J it is necessary 
to include a large number of aggregates (large k) in order to reproduce the absorbance of the 
emulsion for t — 30 s. 

If coalescence occurs, the rate constants determined through the turbidity measurements also 
sample this process. The effect of coalescence should be more pronounced as the salt concen- 
tration increases. In the extreme case in which only attractive forces prevail and coalescence 
occurs instantaneously as the drops collide, the collisions lead to larger drops. The number of 
drops should change as predicted by Smoluchowski for the number of aggregates (Eq. I2D with 
the sole difference that the aggregates of size k are now bigger drops with the equivalent volume 
of k primary drops (Eq. [22). Model M2 represents this case. Curiously, it completely fails to 
reproduce the data corresponding to fast aggregation. 

Model M3 represents the situation in which some collisions between the drops lead to floc- 
culation while others lead to coalescence. The turbidity of an emulsion results from the contri- 
butions of conventional aggregates, bigger drops and mixed aggregates (see Appendix B). If the 
number of mixed aggregates is small, the turbidity can be expressed as a combination of large 
drops and conventional aggregates (Eqs. [28] [29] and 1301 . When the salinity is low and only 
aggregation occurs, g — >1.0. If g is assumed to be equal to 1.0 in Eq. [30] the equation of the 
turbidity for Ml is recovered. Hence, model M3 is capable of reproducing the phenomenon of 
aggregation in the absence of coalescence. As can be inferred from Figs. [T2]and[T3] model M3 
adjusts the turbidity data corresponding to the slow and fast aggregation. Unfortunately, the frac- 
tion of aggregates represented by coefficient "g", surpassed the value of 1.0 for 380 < [NaCl] < 
390 mM. This indicates that the number of mixed aggregates cannot be neglected, and that their 
average effect is being accommodated within the contributions of the conventional aggregates 
and the bigger drops. This can only be possible if the form factor of the mixed aggregates is 
somewhere in-between these two contributions, otherwise the fitting of the experimental data 
will fail if mixed aggregates exist. In this regard, it should be noticed that Eq. [30] can be cast 
into the form of Eq. [31] This last equation is similar to Eq. [34] if the expression inside the 
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brackets is regarded as an average cross section. The cases in which g < 1 correspond to fast 
aggregation (see Table[5]l. In these systems, between 60% and 80% of the primary drops are as- 
sembled into conventional aggregates, and the rest correspond to bigger separate drops resulting 
from coalescence. 

In order to deduce Eq. [30] analytically, it was assumed that the fraction of aggregates (x^ a ) 
and the drops (xk tS ) is the same for each particle size (x^.a = x a , Xk iS - x s ). This requirement is 
rather difficult to attain in practice. On the one hand, it can happen for example, if coalescence 
occurs during binary collisions only, and intra-cluster coalescence does not occur. Hence, only 
a fraction of the collisions will lead to larger aggregates, and the drops will grow by stepwise 
addition of primary particles. On the other hand, mixed doublets are commonly observed in 



macroemulsions 11731. 19711 . This is an indication that either intrac luster coalescence or hetero- 



aggregation occurs for radii of several microns. 

It is also possible that the fraction of aggregates and spheres is similar but not equal at each 
particle size (x^ a ~ x a , Xk, s ~ x s ). In this case, the condition < g < 1 will not be fulfilled, but 
Eq. [30] will still be approximately valid. However, in the general case mixed aggregates cannot 
be disregarded under RLCA conditions. We attempted to solve Eq. [58] for [NaCl]=395 mM, 
using Eqs. [32]or[33]as approximate cross sections. Only the use of Eq. 1331 fitted the data. But 
even in this case, x a < 0, and x m < 0, although x a + x s + x m — 1 .02. Therefore, it is clear that the 
reason why Eq. [58]is successful in reproducing the experimental trend deserves further research. 

Model M4 considers a different aspect of the problem. If the process of coalescence is 
slow enough, the form factor observed by the incoming light of the spectrophotometer might 
correspond to a transient structure, in-between a doublet and a bigger spherical drop. If the 
cross-section of the transient structure is approximated as a sum of spherical and conventional 
aggregate contributions, the expression of the turbidity corresponding to model M3 is recovered. 
Consequently, we tested the analytical form given by Eq. [32] Unfortunately, the results of model 
M4 fitted all experimental data with a similar accuracy than model M3. Hence, the experimental 
data does not allow to discriminate the physical origin of the fittings. Since the use of cr trans 2 (Eq. 
1331 in model M4 (Eq. [34b reproduces the equation of the turbidity of model M3 (Eq. [30b , these 
theoretical approaches rather represent two different interpolations of the cross-section of mixed 
aggregates occurring in emulsions (see Eq. I3TV It can be affirmed that only those models which 
consider the process of partial coalescence can fit the experimental data corresponding to a high 
ionic strength. Moreover, according to the results of the previous paragraph, Eq. [33]appears to be 
a better approximation for the optical cross section of the aggregates than Eq. [32] Despite these 
limitations and the ones previously mentioned regarding the effect of polydispersity, the mixed 
flocculation-coalescence rates observed at high ionic strength (Tables [5] and [6]l are of the same 
order of magnitude as the one predicted by Emulsion Stability Simulations (ESS) for the c ase o f 
very dilute emulsions in the absence of a repulsive force is kf - kfc = 6.4 x 10~ 18 irPs' 1 ll 1 0311 . 
Moreover, the k^ ast obtained are in most cases, half the value of k\\ as predicted by the theory of 
Smoluchowski. 

Fig. [16] shows a comparison of all stability ratios calculated in this work. The ones corre- 
sponding to models M3 and M4 were calculated using the flocculation rates of Tables [5] and [6] 
In our opinion the agreement is satisfactory. 

Finally, it is necessary to discuss the effect of polydispersity, given our emulsions are rather 
polydisperse ( CV - 35%). The influence o f poly dispersity has been investigated by many autho rs 
(Mullholland EH, Friedlander and Wang [fl47l[l48h . Cohen and Vaughan ThI, and Lee lfl50ll ). 
The continuous distribution of particle sizes is a non-linear integro-differential equation whose 
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kernel depends on the size of the colliding aggregates. Meesters and Ernst [151] demonstrated 
that the use of a constant kernel in this equation leads to an exponential frequency distribution 
of particle sizes. Smoluchowski result approaches this exponential distribution at long times. If 
a homogeneous kernel is used, and the degree of homogeneity is greater than zero, bell-shaped 
distributions are obtained Il20ll . In this case, collisions between large and small aggregates occur 
more rapidly than in the case of a constant kernel. This extra depletion of smaller aggregates 
leads to bell-shaped particle size distributions. 



In regard to polydispersity the work of Lee 111501 is especially relevant to this investigation 



This author assumed that the initial and subsequent drop-size distributions could be represented 
by a log-normal function. In this case, an expression for the number of aggregates produced by 
flocculation -similar to the one of Smoluchowski- can be deduced: 

(37) 



(1 + CpkfHot) 
where: 

C p = 1 + expO 2 ^)), (38) 

and erg = Tr(t — 0). 

Coefficient C p corrects the flocculation kernel for polydispersity. According to Lee Il50ll the 
size distribution approaches that having a geometric standard deviation of (fa, = 1.32 regardless 
of the value of cF . The difference between the predictions of Eq. [2] and the ones of Eq. [37] 
for an initially monodisperse system is less than 4%: C p = 1.08 115011 (C p =1.13 according to 



Friedlander and Wang [147]). However, if the initial distribution is polydisperse, the difference 
between both equations becomes much larger, and increases with the value of ao- Even in this 
case the asymptotic limit is approached, but at considerably longer times. In any event, it should 
be noticed that coefficient C p multiplies kf. Hence, a fit of the data to Eq. [2]leads to an uncertainty 
in the value of kj which increases with the initial polydispersity of the system. 

It is clear from Figs. [17] and [18] that the increase of the polydispersity of the initial (mother) 
emulsions decreases the quality of the fittings. However, model M3 is less susceptible to this 
variable than model Ml. This is probably due to its additional fitting parameter (g). In any 
case, the introduction of polydispersity does not change the curvature of the absorbance curve, 
supporting our main results that models Ml and M2 do not describe the process correctly at high 
salinity, whereas M3 and M4 fit the data correctly. 



6. Conclusions 

The stability of a dodecane/water nanoemulsion stabilized with SDS was evaluated using 
turbidity measurements. For this purpose, the variation of the turbidity as a function of time 
was expressed in terms of the flocculation rate of the dispersion. Four approximate theoretical 
models were postulated in order to represent the contribution of the aggregates of the emulsion 
to the turbidity. The complete set of experimental data was only reproduced by those models 
which take into account the existence of the conventional aggregates as well as bigger spherical 
drops resulting from coalescence. The adjustment of the experimental turbidity to the corre- 
sponding theoretical expressions allowed the determining of the value of the mixed flocculation- 
coalescence rate. The stability ratios (W) evaluated using these rates are consistent with those 
determined experimentally from the initial slope of the absorbance as a function of time. 
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As previously pointed out, the complexity of the emulsions results from the coupling of 
several destabilization mechanisms. In this report, a "window" of compositions where the pro- 
cesses of flocculation and coalescence outweigh significantly other destabilization mechanisms 
is found. It is important to remark that the mixed flocculation-coalescence rates shown in Tables 
|5]and|6]do not correspond to the original nanoemulsion synthesized. They only appraise the ef- 
fect of the referred processes in the presence of a considerable amount of salt. Several technical 
difficulties appear when the salt concentration is lowered much further (< 300 mM). Hence, it is 
not generally possible to evaluate the flocculation rate of the original emulsion. This restriction 
also applies in the case of suspensions. 
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8. Appendices 

8.1. Appendix A: Alternative procedure for evaluation ofEqs. \T4\\T5\\T8\ and \19\ 

It is clear from Eq. [18]that the integr als o f the right hand side do not depend on k. Thus, the 



total cross-section, o~k can be written as [133]: 



a k = ko-y + -na 2 a\m - l) 2 ^ 7 y , (39) 



Iij = I n = f { f "... f P(n, ...^. /'"^f 71, "" y " ) dyi...dy„\ x W)(l + cos 2 V)sin&d§, 
Jo U J hs u (y u ...,y n ) I 

(40) 

where p{y\ , y„) is a flat distribution function, and n is the number of particles between particles 
; and j. 



In the case that we consider a single particle (n = — 1), Eq. |40]can be written as 113311 : 



Jo 



/_! = | Pi(#)(l +cos 2 ff)sinM§. (41) 



For the doublet, particles are directly connected and n = 0, simplifying Eq. l40lto: 

7 = 2 f Sm( ~ hS ' 12) TMffXl + cos 2 &)sin$d&. (42) 

Jo '"12 

For the case n — 1, we have: 



r r 2 *' 3 sin(/zs 13 (ri)) , 

7i=2 p , \ P l (&)(l+cos 2 $)sin§d&. (43) 

Jo J-2;r/3 hsnijx) 
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Finally for n = 2, we can write: 

-2x/3 J-2?r/3 



/ 2 = 2 f f*V^^P.(W + cos^sinMO. (44) 

Jo J-2;r/3 J-2;r/3 «.Sl4(yi>72) 



Using the definition of /„, the total cross-section of an aggregate of N particles can be written 
as follows: 

cr„ = No-i + -na 2 a\m - l) 2 J^(N -n - l)/„, (45) 

since there are N — n - 1 pairs of particles with n particles between them. Here, kcs I-\ — o~\ and 
k cs = (4/9)7TflV(m - l) 2 . 

Now, using Smoluchowski kinetics Q\ and considering a unique kinetic constant for coagu- 
lation, kf independent of the aggregate form or size, the change in the number of aggregates with 
k particles is given by: 



drik(f) 1 



dt = X n <'W"/W _ k f n ^ Zi ni(f) - (46) 

•J=l, i+i=k i'=l 

Thus, fe-particle aggregate concentration as a function of time, n*(f) can be obtained analytically 

HQ 

t k-\ 

(^ T *p) 

where «o is the initial concentration of particles. The turbidity can be expressed as: 

Puertas and de las Nieves lfl3l demonstrated that this expression retains the most important 
contributions to the summation if it is truncated at n — 2. 

8.2. Appendix B: Balance of aggregates for a flocculating-coalescing emulsion 
In the case where only aggregation occurs, the balance of mass gives: 



^ n* = n tota i. (49) 

k 



The particles can increase their size due to aggregation, coalescence, and aggregation of 
bigger drops with primary particles to form "mixed" aggregates Fig. [19] 

ntotai = »i + 2_u + + ' ( 5 °) 

k=2 

where, n^s is the number of the drop resulting from coalescence, n^ a represents the number of 
aggregates with k primary particles, and n^ m is the number of mixed aggregates resulting from 
the coalescence and aggregation of k particles. 
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If both aggregates and spherical drops form as a result of binary collisions, they should 
be a fraction (x) of the total number of aggregates predicted by Smoluchowski for the case of 
irreversible flocculation: 

nk, s = n k Xk, s , (51) 
rik, a = n k x kta , (52) 
nk,m = n k Xk, m . (53) 



)■ (54) 



jfc=2 i=2 



Therefore, 
For each £: 

^ + JCi, a + % m = 1. (55) 

Hence, 

That is, although the number of initial particles is different from the number of final particles, the 
total number of aggregates is the same for a flocculating suspension and a fiocculating-coalescing 
emulsion. These two systems differ in the types of aggregates formed. When coalescence oc- 
curs, it generates mixed aggregates which are counted as aggregates of size k in the absence of 
coalescence (their total volume corresponds to the sum of k primary particles). Therefore, the 
number of aggregates does not change with respect to the definition of Smoluchowski, but their 
structure does. 

Now, the turbidity can be redefined as: 

t — ii\ o~\ + V x Ks n k cr k ^ + V x kta n k a-k,a + V x kim n k cr^ m , (57) 

k=2 k=2 k=2 

t = Hi<xi + ^ n k (x k , s cr ktS + x Ka <Tk,a + Xk, m Ck,m) • (58) 



k=2 

If n k ^ m is negligible (x k , m = 0), then 

T = nicr x +^n k (x kiS cr Ks + x Ka cr kA ) . (59) 

k=2 

In this case, Xk, s + x k ^ a = 1, therefore Eq. |59]is equal to: 

T = n\<T\ + ^ Uk [(1 - X k , a ) CT k , s + X Ka O-^a\ . (60) 
k=2 

If all x kM — x a and x k , s — x s , then 

t = n\(T\ + 2^ n k [(1 - x a )o- k , s + ]■ (61) 

k=2 
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If it is found that Eq. |6T]is confirmed by the experiments but x s + x a + 1, this means that the 
average contribution of mixed aggregates is being incorporated into cr^ a and cr^ s . If we equalize 
Eq. [60]to the exact form of the turbidity (Eq. [58}, r = r = T exp , then 

/ \ nk (Xk, s O~k,s + Xk,aCTk,a + Xk, m O~k,m) = ^ «t )■ (62) 

Considering Eq. [55] we have: 

It is clear then that the fact that x s + x a + 1 comes from mixed aggregates contributions. 
Moreover, this contribution can be evaluated since x, s and x, are known. For this purpose, it is 
necessary to provide an approximate expression for cr^m- 

8.3. Appendix C: Calculation of errors 

In order to assign an error to the theoretical estimations of kf from models Ml to M4, an av- 
erage error of the turbidity was approximated by the average difference between the experimental 
values of the turbidity (j exp ) and the theoretical prediction (t/, ( ). 

The average error of the theoretical kinetic constant (kf) was calculated from the deviates of 
Abs (AAbs) 

AAbs = — — V \Abs fh - Abs exp \ , (64) 

Ntotal 

Ntotai represents the number of values of Abs for each experimental turbidity curve. The error of 
kf (Akf) was calculated from the propagation of errors of theoretical form of r for each model, 
obtained using symbolic algebra (Mathematica 8.0.1.0): 

AAbs 

M ' = WIti7\ 

Error bars of the initial slopes showed in Fig. |6]correspond to the standard deviation of three 
measurements. With these errors and the propagation of errors of Eq. [9] A£n was obtained. 
The error of critical coagulation concentration (ACCC) is given by: 

(66) 

The deviation of slope (Ac) and deviation of the intercept (Ab) were determined from the least- 
squares fitting of the W vs. [NaCl] curves. 

8.4. Appendix D: Comparison between Eq. \14\and Eq. \39\ ( Model Ml ) 

A comparison between the two forms of evaluation of Ml is shown in Fig. [20] In the case 
of Eq. [14] a value of k max = 100 was used. As can be seen, the adjustment of Eq. [14] to the 
experimental data is very good at the initial stage. However, the final points (long times) do not 
fit the data. On the contrary, Eq. [39]fits very well in the final stage, while the initial points escape 
the fitting. 

It is necessary to point out that the quality of the fitting depends on the values of the average 
radius of primary particles. It is well known that distinct average radii can be obtained using 
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size distributions weighted in intensity, volume, area or number of particles. In the present 
calculations, we chose an "optimal" radius (R/i,) within the size range R exp + AR exp (Table |2]l 
previously measured, in order to maximize the agreement between the experimental data and the 
theoretical models M3 and M4. Thus, in order to compare all the models on a similar basis, the 
same set of radii (/?/,>) was considered to evaluate model Ml using Eqs. [14]and[39] Model Ml 
does not fit the DLCA data. The values of k f and f for 380 < [NaCl] < 395mM predicted from 
both equations are shown in Table [8] It is clear that the results of Eq. [14] are of lower quality. 
Therefore, Eq. |39]was used for calculating the kinetics constants of model Ml (Table|4| under 
DLCA conditions. 
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Table 1 : Composition of the systems employed. 



Parameters 


Starting 


Concentrated 


Final 




system 


nanoemulsion 


nanoemulsion 


<P 


0.84 


6.5 X 10" 2 


3.2 X 10" 4 


[SDS] (M) 


0.35 


2.4 X 10" 2 


8.0 X 10" 3 


[NaCl\ (M) 


1.40 


2.0 X 10" 2 


9.5 x 10" 5 


[Isopentanol] (wt%) 


6.50 


4.6 x 10-' 


2.2 x 10" 3 


Droplets/mL 




4.1 x 10 13 


2.0 x 10" 



Table 2: Characterization of the nanoemulsions. 



Mean radius (nm) 72.5 
Standard deviation (nm) 25.4 
Coefficient of variation (CV) 35 % 
I -potential (mV) -81.7 ± 2.5 



Table 3: Values of fen from Eq. [9] 



[NaCl] R exp k n 
(mM) (nm) (mV) 



380 


75.2 


(2.1 ±0.8) X 10" 


■19 


390 


78.0 


(3.5 ± 0.6) X 10" 


•19 


395 


74.3 


(3.8 ± 1.9) x 10" 


19 


400 


75.1 


(1.2 ± 0.1) x 10" 


■18 


410 


78.0 


(1.6 ± 0.6) x 10" 


■18 


420 


76.1 


(1.9 ± 0.2) x 10- 


■18 


430 


74.3 


(3.0 ± 0.4) x 10- 


■18 


450 


74.3 


(2.8 ± 0.3) x 10" 


■18 


475 


75.0 


(3.2 ± 0.6) x 10- 


■18 


500 


72.3 


(3.3 ± 0.3) x 10- 


18 


550 


72.3 


(3.5 ± 0.7) x 10- 


■18 


600 


75.0 


(3.1 ±0.1) x 10- 


■18 
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Table 4: Kinetic constants obtained from model Ml. 



[NaCl] 


Rfit 




tf),exp 


tojit 


(mM) 


(nm) 


(m 3 .?" 1 ) 


(s) 


(s) 


380 


64.0 


(1.4 ± 0.3) X 10" 19 


11.9 


11.0 


390 


64.0 


(3.3 ± 0.2) X 10" 19 


7.2 


7.0 


395 


65.0 


(3.3 ± 2.5) x 10" 19 


6.3 


7.8 


400 


66.0 


(4± 18) X 10~ 18 


6.3 


7.5 


410 


66.0 


(4±38)x Kr 18 


6.2 


6.5 


420 


75.0 


(2±44)x 10~ 17 


7.4 


7.5 


430 


63.0 


(6±70)x 10~ 18 


6.3 


7.0 


450 


69.3 


(5 ± 27) X 10~ 18 


6.5 


6.7 


475 


58.0 


(5±35)x 10~ 18 


7.3 


7.4 


500 


69.7 


(7±2)x 10~ 18 


6.4 


6.6 


550 


72.7 


(7 ± 29) x lO" 19 


6.4 


8.9 


600 


57.7 


(5 ± 148) x 10" 19 


7.3 


7.8 



Table 5: Kinetic constants obtained from model M3. 



[NaCl] 


R fit 


k/H3 




g 


k),exp 


tQ, fit 


(mM) 


(nm) 


(mV 1 ) 






(s) 


(s) 


380 


64.0 


(1.7 ± 0.1) x 10" 


•19 


1.35 


11.9 


13.0 


390 


64.0 


(2.2 ± 0.1) x 10" 


-19 


1.07 


7.2 


6.6 


395 


65.0 


(3.5 ± 0.3) X 10" 


■19 


1.00 


6.3 


7.1 


400 


66.0 


(6.7 ± 0.2) X 10" 


■19 


0.62 


6.3 


7.2 


410 


66.0 


(7.9 ± 0.4) X 10" 


■19 


0.64 


6.2 


6.8 


420 


75.0 


(2.2 ± 0.1) X 10- 


■18 


0.69 


7.4 


7.6 


430 


63.0 


(1.0 ± 0.1) X 10- 


■18 


0.78 


6.3 


7.1 


450 


69.3 


(1.3 ± 0.1) x 10- 


-18 


0.73 


6.5 


6.9 


475 


58.0 


(8.4 ± 0.2) X 10- 


-19 


0.79 


7.3 


7.5 


500 


69.7 


(1.5 ± 0.1) x 10- 


-18 


0.78 


6.4 


6.8 


550 


72.7 


(1.5 ± 0.1) x 10- 


-18 


0.81 


6.4 


8.9 


600 


57.7 


(6.7 ± 0.1) x 10- 


-19 


0.81 


7.3 


7.4 



Table 6: Kinetic constants obtained from model M4. 



[NaCl] 


Rfit 






P 




tojit 


(mM) 


(nm) 


(mV) 






(s) 


(s) 


380 


64.0 


(2.0 ± 0.1) X 10- 


-W 


1.68 


11.9 


13.1 


390 


64.0 


(2.2 ± 0.1) X 10" 


-19 


1.11 


7.2 


6.6 


395 


65.0 


(3.5 ± 0.3) X 10- 


-19 


1.00 


6.3 


7.1 


400 


66.0 


(6.1 ±0.1) x 10- 


-19 


0.43 


6.3 


7.2 


410 


66.0 


(7.0 ± 0.3) X 10- 


-19 


0.45 


6.2 


6.8 


420 


75.0 


(1.9 ± 0.1) X 10- 


-18 


0.47 


7.4 


7.5 


430 


63.0 


(9.0 ± 0.2) X 10- 


-19 


0.52 


6.3 


6.8 


450 


69.3 


(1.2 ± 0.1) x 10- 


-18 


0.55 


6.5 


6.9 


475 


58.0 


(7.3 ± 0.1) x 10- 


-19 


0.61 


7.3 


7.4 


500 


69.7 


(1.4 ± 0.1) x 10- 


-18 


0.60 


6.4 


6.9 


550 


72.7 


(1.4 ± 0.1) x 10- 


-18 


0.67 


6.4 


8.9 


600 


57.7 


(5.9 ± 0.1) x 10- 


-19 


0.65 


7.3 


7.3 
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Table 7: Values of critical coagulation concentration. 



Data 


CCC (mM NaCl) 


From W exp 


412 ±6 


From Wjt n 


411 ± 3 


From Wm3 


415 ± 9 


From Wm4 


418 ± 13 



Table 8: Comparison between the values of k/ and tg for 380 < [NaCl] < 395 mM using Eq. l39l and Eg. 1 141 



[NaCl] 


k f -Eq.M 


tQJit 


k f -Eq.[39\ 


tajtt 




(mM) 


(mV 1 ) 


(s) 


(mV 1 ) 


(s) 


(s) 


380 


7.52 X 10~ 2U 


8.3 


1.37 X 10~" 


10.9 


11.9 


390 


1.93 x 10~ 19 


6.4 


3.32 X lO" 19 


7.0 


7.2 


395 


3.82 X 10~ 19 


6.6 


3.30 X lO" 19 


7.8 


6.3 



Figure Captions 

Figure 1: Definition of y for a chain of three particles (a) and four particles (b). 

Figure 2: [NaCl] vs. [SDS] phase diagram of the water/NaCl/isopentanol/SDS/dodecane system 
at 6.5 %wt isopentanol and f w = 0.20. 

Figure 3: Droplet size distributions of three independent nanoemulsions synthesized from three 
different mother emulsions measured with a Malvern Zetasizer 2000. 

Figure 4: Evolution of R exp as a function of time for several nanoemulsions. 

Figure 5: Typical curve of Abs vs. t. 

Figure 6: Aggregation kinetics of the nanoemulsion. 

Figure 7: Results of two measurements of the initial droplet size distribution of the same na- 
noemulsion using a BI-200SM Goniometer: (a) Contin, (b) Cumulants analysis. 

Figure 8: Temporal evolution of the droplet size distributions for a nanoemulsion with 380 mM 
NaCl using a BI-200SM Goniometer: (a) Contin, (b) Cumulants analysis. 

Figure 9: Temporal evolution of the droplet size distributions for a nanoemulsion with 600 mM 
NaCl using a BI-200SM Goniometer: (a) Contin, (b) Cumulants analysis. 

Figure 10: Droplet size distributions after 1 min of addition of different concentrations of salt 
using a BI-200SM Goniometer: (a) Contin, (b) Cumulants analysis. 

Figure 11: Comparison of the theoretical curves of absorbance predicted by models Ml and M2 
with the experimental measurements of slow aggregation ([NaCl] = 390 mM). 

Figure 12: Comparison of the theoretical curves of absorbance predicted by models M3 and M4 
with the experimental measurements of slow aggregation ([NaCl] = 390 mM). 

Figure 13: Comparison of theory and experiment for the case of fast aggregation ([NaCl] = 600 
mM). 
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Figure 14: Predictions of model Ml for k max = 2, 3,4 and 100 ([NaCl] = 390 mM). 
Figure 15: Comparison between W*,, andW exp . 

Figure 16: Comparison between theoretical and experimental values of W. 

Figure 17: Comparison of the behavior of T aprox with t(R fu) for models Ml and M3 at [NaCl] = 
390 mM. (The number of data points of the experimental curve was reduced in order to 
highlight the quality of the fittings). 

Figure 18: Comparison of the behavior of T aprox with t(R fn) for models Ml and M3 at [NaCl] = 
475 mM. (The number of data points of the experimental curve was reduced in order to 
highlight the quality of the fittings). 

Figure 19: Types of aggregates formed in an emulsion. 

Figure 20: Evaluation of the model Ml using Eq. [14] and Eq. [39] 
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Figure 1: Definition of y for a chain of three particles (a) and four particles (b). 



- w + o 

t m | 

— a i 

A , 


A 
A 


o a 

m ^ ^ 
- - ' A 


A 


— W+D+O ! 

- e ! 




W+LC+O 




— ' 
v 1 






* 




» A 
A- 


A 


A 


w+o 

m 




A 


A 




A 


Equilibrated systems 
Starting system 




" 


1 





1 1 



3456789 10 11 

[SDS] (wt%) 



Figure 2: [NaCl] vs. [SDS] phase diagram of the water/NaCl/isopentanol/SDS/dodecane system at 6.5 %wt isopentanol 
and f w = 0.20. 
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Figure 3: Droplet size distributions of three independent nanoemulsions synthesized from three different mother emul- 
sions measured with a Malvern Zetasizer 2000. 
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Figure 4: Evolution of R exp as a function of time for several nanoemulsions. 
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Figure 5: Typical curve of Abs vs. t. 
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Figure 6: Aggregation kinetics of the nanoemulsion. 
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Figure 7: Results of two measurements of the initial droplet size distribution of the same nanoemulsion using aBI-200SM 
Goniometer: (a) Contin, (b) Cumulants analysis. 




R (nm) R (nm) 

Figure 8: Temporal evolution of the droplet size distributions for a nanoemulsion with 380 mM NaCl using a BI-200SM 
Goniometer: (a) Contin, (b) Cumulants analysis. 
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Figure 9: Temporal evolution of the droplet size distributions for a nanoemulsion with 600 mM NaCl using a BI-200SM 
Goniometer: (a) Contin, (b) Cumulants analysis. 




Figure 10: Droplet size distributions after 1 min of addition of different concentrations of salt using a BI-200SM Go- 
niometer: (a) Contin, (b) Cumulants analysis. 
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Figure 11: Comparison of the theoretical curves of absorbance predicted by models Ml and M2 with the experimental 
measurements of slow aggregation {[NaCl] = 390 mM). 
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Figure 12: Comparison of the theoretical curves of absorbance predicted by models M3 and M4 with the experimental 
measurements of slow aggregation ([NaCl] = 390 mM). 
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Figure 15: Comparison between Wt n and W el 
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Figure 16: Comparison between theoretical and experimental values of W. 
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Figure 17: Comparison of the behavior of T aprax with t(R/, ( ) for models Ml and M3 at [NaCl] = 390 mM. (The number 
of data points of the experimental curve was reduced in order to highlight the quality of the fittings). 




Figure 18: Comparison of the behavior of r aprox with T(Rfu) for models Ml and M3 at [NaCl] = 475 mM. (The number 
of data points of the experimental curve was reduced in order to highlight the quality of the fittings). 
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Mixed aggregate (m) 



Figure 19: Types of aggregates formed in an emulsion. 



47 




48 



